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Abstract 

A general system of particles (of one or several species) on a one dimen- 
sional lattice with boundaries is considered. Two general behaviors of such 
systems are investigated. The stationary behavior of the system, and the 
dominant way of the relaxation of the system toward its stationary sys- 
tem. Bases on the first behavior, static phase transitions (discontinuous 
changes in the stationary profiles of the system) are studied. Based on the 
second behavior, dynamical phase transitions (discontinuous changes in 
the relaxation-times of the system) are studied. The investigation is spe- 
cialized on systems in which the evolution equation of one-point functions 
are closed (the autonomous systems). 
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1 Introduction 



The study of the reaction-diffusion systems, has an attractive area. A reaction- 
diffusion system consists of a collection of particles (of one or several species) 
moving and interacting with each other with specific probabilities (or rates in the 
case of continuous time variable). In the so called exclusion processes, any site 
of the lattice the particles move on, is either vacant or occupied by one particle. 
The aim of studying such systems, is of course to calculate the time evolution 
of such systems. But to find the complete time evolution of a reaction-diffusion 
system, is generally a very difficult (if not impossible) task. 

Reaction-diffusion system have been studied using various methods: ana- 
lytical techniques, approximation methods, and simulation. The success of the 
approximation methods, may be different in different dimensions, as for exam- 
ple the mean field techniques, working good for high dimensions, generally do 
not give correct results for low dimensional systems. A large fraction of analyt- 
ical studies, belong to low-dimensional (specially one-dimensional) systems, as 
solving low-dimensional systems should in principle be easier ,lli2„3||4i,5|i6j|71|Hl 

EmniiiiMTiini. 

Various classes of reaction-diffusion systems are called exactly-solvable, in 
different senses. In jJJ] and J^, integrability means that the A^-particle con- 
ditional probabilities' S-matrix is factorized into a product of 2-particle S- 
matrices. This is related to the fact that for systems solvable in this sense, there 
are a large number of conserved quantities. In (16 , 17, 18 , 19 , 20.21 .22..23,..24il25| . 
solvability means closedness of the evolution equation of the empty intervals (or 
their generalization). 

Consider a reaction-diffusion system (on a lattice) with open boundaries. 
By open boundaries, it is meant that in addition to the reactions in the bulk of 
the lattice, particles at the boundaries do react with some external source. A 
question is to find the possible phase transitions of the system. By phase tran- 
sition, it is meant a discontinuity in some behavior of the system with respect 
to its parameters. Such discontinuities, may arise in two general categories: in 
the stationary (large time) profiles of the system, and in the time constants de- 
termining the evolution os the system. In the first case, static phase transitions 
are dealt with; in the second case, dynamical phase transitions. As mentioned 
before, the task of finding the complete evolution of a general reaction-diffusion 
system is very difficult (if not impossible) . So our studies are limited to a certain 
class of reaction-diffusion systems, and in them to certain properties. To be spe- 
cific, we deal with one-dimensional systems for which the evolution equation of 
one-point functions (probabilities of finding a particle of a certain kind on a cer- 
tain point) is closed [autonomous systems). For these systems, we only consider 
the final (stationary) profile of the one-point function, and the relaxation-time 
(the most significant one) towards the stationary profile. 

In a ten-parameter family of one-species reaction-diffusion processes 
with nearest-neighbor interaction was introduced, for which the evolution equa- 
tion of n-point functions contains only n- or less- point functions. The aver- 
age particle-number in each site has been obtained exactly for these models. 
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In [271|5S], this has been generahzed to multi-species systems and more-than- 
two-site interactions. In IJHEOllSElli the phase structure of some classes 
of single- or multiple-species reaction-diffusion systems have been investigated. 
These investigations were bases on the one-point functions of the systems. 

The scheme of the present article is as follows. In section 2, the condi- 
tions for the most general multi-species reaction-diffusion models with nearest- 
neighbor interactions, to be solvable, or autonomous, are obtained. In section 3, 
autonomous single-species reaction-diffusion models with boundaries have been 
investigated. It is shown that changing the bulk rates may lead to a static phase 
transition, while changing the bulk or boundary rates may lead to a dynamical 
phase transition. In section 4, as an example, an asymmetric generalization of 
the Glauber model at zero temperature is introduced. This system exhibits both 
static and dynamical phase transitions. In section 5, as an example of systems 
with more than two-site interaction, Glauber model at finite temperature has 
been considered. With specific boundary interactions, the Glauber model with 
open boundaries is also autonomous. It is shown that although this system does 
not show a static phase transition, it does exhibit a dynamical phase transition. 
In the last section, as an example of multi-species reaction-diffusion systems, an 
asymmetric generalization of the zero temperature Potts model is introduced. 
This system exhibits both a static and a dynamical phase transition. 

2 Autonomous reaction-diffusion systems, with 
nearest-neighbor interactions 

In |27|. multi-species reaction diffusion with nearest-neighbor interactions on a 
one-dimensional lattice are studied, and criteria are obtained that such systems 
be autonomous, meaning that the evolution equation for the one-point functions 
be closed. The summary of the work follows. 
Let the Hamiltonian of the system be 

L-l 

H=^ff,..+i, (1) 
1=1 

where 

:= 1(g) ■ ■ ■(g) l (g)g(g) 1(g) ■ ■ ■(g) 1 . (2) 

■i-l L-i-l 

1 is the {p+ 1) X (p+l) unit matrix, and H is a. {p + 1) x (p + 1) matrix (for a 
p species system) , the nondiagonal elements of which are nonnegative. H also 
satisfies 

sH^O, (3) 
where s is a covector with all of its components equal to 1: 

Sa = 1. (4) 
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The number operator of a particle of type Aa in the site i is denoted by nf, 
where Ap+i can be the vacancy. These number operators are of the form 

1 (g) • • • (g) (g) 1 (g) • • • (g) 1, (5) 

i — 1 L—i 

where n" is a diagonal (p+1) x matrix, the only nonzero element of which 

is the element a of the diagonal. It is clear that the vector n whose components 
are n", satisfies 

sn=l, or San°' — i- (6) 

The Hamiltonian Ti. is the generator of time translation, by which it is meant 
that 

^P{t)=HPit\ (7) 

where P is a the vector in the (p+1)^ dimensional space, the component p^i - ^^- 
of which is the probability of finding a particle of type Aa^ in site 1, etc. Clearly, 
the components of this vector are nonnegative and satisfy 

SP-1, (8) 

where 

S := s (g • • • (g s . (9) 



L 

The fact that P should have these properties for all times, means that the 
nondiagonal elements of H should be nonnegative and that the action of Tl 
(from right) on S should be zero. But the nonnegativity of the nondiagonal 
elements of H, and the condition , are sufficient for these two conditions on 
H to be satisfied. 

The expectation value of any operator O is 

{0) = SOP. (10) 

This is also true for nf, from which it is seen that 



-«) = S<HP. (11) 

Now the question is, under what conditions the right-hand side of ((TT|l can be 
written in terms of the one-point functions n^^sl 

To answer this, one notes that nf commutes with all of terms in Ti, except 
possibly with and Hi^i^i. Using 

s g) s (an) g) (bn) H = (a g) b i/)Q;3 s g) s n" g) n'', (12) 

which is true for any two covectors a and b, and using (jSJ, it is seen that 

^«) = <) + s, H^\, « <+i). (13) 
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The right hand side is expressed in terms of two-point functions not one-point 
functions. However, from © it is seen that the two-point functions are not 
independent. In order that an expression B^i, be expressible in terms of 

the first power of number operators, it is necessary and sufficient that 

Bf,^ = iB^^s^, + Sf,2B^. (14) 

So the necessary and sufficient condition for the system to be autonomous, is 
that 

+ (15) 

where 

Irra — rra7 

• — '='7 ^U: 

(16) 

Note that even if pSI) is satisfied, jH^s are not uniquely determined through 
which. One can change them like 

}i?"^^}i?"^+>"s;„ (17) 

with 

E>-0- (18) 

3 

For the simplest case p ~ 1 (single-species systems) (|15|l consists of two inde- 
pendent constraints, leaving a ten-parameter family of autonomous systems out 
of the twelve-parameter family of the general systems. 

The condition (|15|l for a system to be autonomous, can be easily extended to 
systems with interaction-ranges longer than the nearest neighbor. This has been 
done in [^Hj. The result is that if the interaction is in a block of k neighboring 
sites, then one should have 

'H^p,...p, - \H^p, Sft • • • + • • • + +sp, ■ ■ ■ sp,_, IH^0,, (19) 

where 

:= • • • S7,_, s^,^, • • • s^, H'^--"- (20) 

Assuming that (|15|l holds, one can write the evolution equation for the one- 
point functions as 

^{n,)=lH{n,_^) + {lH + \H){n,)+\H{n,+^), l<i<L. (21) 
at 

For i = 1, at the right-hand only the last two terms remain; for i — L, only the 
first two terms. This is a linear differential-difference equation for the vectors 
(n). 
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3 Autonomous single-species reaction-diffusion 
systems with boundaries 

For a single-species system, the condition (|15|) can be rewritten in the more 
explicit form 

rrOl rrOO , rrOl , rrOO rrll rrlO , ttU , rrlO n 

— li 11 — n 11 + U 10 + U 10 — U 01 — ^ 01 + -ti 00 + 00 — • t), 

rrlO TirOO rrll titOI , rrlO , titOO , rrll , rrOl n 

— li 11 — U 11 — -ti 10 — H 10 + ^ 01 + ^ 01 + ^ 00 + 00 — • U. 

(22) 

Here the state Ai is an occupied site, while is a vacancy. Defining 



u 


ttW 

■ — n 01- 


t -ti 01 


V 


rrOl 
■ — n 10 - 


^H°\o 


u 


■.= H^\i- 


FH""ii 


V 


■■=H^\i- 


Fi/°°ii 


w 


■ ~ H^^oo - 


t ii 00 


s 


H^^oo - 


H 00 


w 


■ = H^^oi - 


H 01 


s 


■— H^^io - 


^ 10 



one can write (|22(l as 

U + S = fi + s 

V + w ^ V + w. (24) 

The evolution equation of (rii), for 1 < « < L is of the form 12111 . For i — 1 
and i = L, only half of the terms at the right-hand side of (|21|l are there. 
Moreover, if there is injection and extraction of particles at the end sites, terms 
corresponding to these rates should also be added. Finally, for the single species 
case, denoting nj by n^, one can write instead of the vector equation 1)21(1 . a 
scalar equation for rii, as n] + = 1. denoting the injection and extraction 
rates at the first site by a and a' respectively, and those at the last site by b and 
6', the evolution equation of the one-point functions can be seen to be 

{hi) = — {v + w + U + s)(n,;) + {v — v){ni+i) + [u — u){ni^i) 

+ w + s, 1 <i < L 
{hi) — — [v + w){ni) + [v ~ v){n2) +w + a{l - (rii)) - a'(ni) 
("L> ^-{u + s){nL) + {u- u){nL-i) +s + h{l- {ul)) - b'{nL), (25) 

3.1 The static phase transition 

The stationary-state to 1251 is 

{rij) = C + Dizl+D2zl (26) 
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where z^'s satisfy 

-r +pzi + qz^^ = 0, (27) 

and we take Z2 to the root of the above equation with larger absolute value, and 
the new parameters p, q, and r are defined through 

p -.—v — V 

q :=u — u 

r :=u -\-S'\'V + w = u-\-s~\-v + w 

=t2 + s + u + ui = + s + -D + iD, (28) 

Noting that the rates (nondiagonal elements of H) arc nonnegative, it is seen 
that 

r>\p + ql\p-q\. (29) 

More over, it is seen that if r = 0, then all the parameters introduced in H23|l 
vanish, and (jij) would be constant (for 1 < i < L). So apart from this trivial 
case, r is positive and in fact can be rescaled to one (by a proper redefinition of 
time). Hence, as long as only the stationary profile of the one-point function is 
considered, there are only two-parameters in the system: 

a; := -, and y := -. (30) 

r r 

The physical region corresponding to these parameters is a square the bound- 
aries of which are |a; ± y| = 1. 

It can be seen that inside the physical square, both of the roots of (|27|l are 
real, one (zi) between —1 and 1, the other (22) outside that region. So in the 
thermodynamic limit L — > 00, 

{n^) + Diz{, j<^L 

{fij) ^ C + D2 zi^^-\ L-j<i. (31) 

C, Di^s, and z^'s are smooth functions of rates, and there is no phase transition. 

However, if |a; -I- y| = 1, then one of the roots of l|77|l would he x + y; that 
is, the absolute value of one of the roots will be 1. The other roots would be 
y/x. If \y/x\ < 1, then zi — 1 and from (|31|l it is seen that near the end of the 
lattice, \{nj) — C| is essentially constant. If \y/x\ > 1, then Z2 — 1 and from 
H31fl it is seen that near the beginning of the lattice, \{nj) — C\ is essentially 
constant. So there is a phase transition at x = y — ±1/2. This static phase 
transition manifests itself as a change in the slope of the profile of the one-point 
function, near the end or the beginning of the lattice. 

3.2 The dynamical phase transition 

The dynamical phase transition is related to the relaxation time of the system 
evolving towards its stationary configuration. To find this relaxation time, one 
writes the homogeneous part of l|25|l as 

N.^h^^Nj, (32) 
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where iVj is the (n^) minus its stationary value. Then, one seeks the eigenvalues 
of h. The eigenvalue with the largest real part, determines the relaxation time: 



1 



Re(£;„ 



(33) 



The eigenvector equations read 



-{v + w + u + s)X.j + (u - v)Xj+i + {u — u)Xj^i, j 7^ 1, L 



EX, 

EXi = -{v + w + a + a')Xi + {v - v)X2, 
EXL^-{u + s + b + b')XL + {u- u)Xl-i. 



A solution for this is 
with 
and 



E 



Xj — a z{ + Pz2, 
{v + w + u + s) + {v — v)zi + [u- u)z^^, 



{v - v){azf + I3zl) -{E + a + a' + v + w){azi + /?Z2) = 
{u - u)(azf ~i + Pz^-^) -{E + b + b' + u + s){az^ + f3z^) = 0. 



(34) 

(35) 
(36) 

(37) 



One puts E from (|36|l in (|36|l . and demands l|36|l to have nonzero solutions for 
a and /3. The result is 



[(u — u) + zi6a\[(v — v)z2'^^ + zl^Sb] — [{u — u) + Z25a\ 



[{v - w)zf + 56] = 0, (38) 



where 



(5a := a + a' — (u + s), 
(56 := 6 + 6' - (w + w). 



Defining the new variables Zi as 



Zi : — Zi 



u — u 



(39) 



(40) 



it is seen that 

E = -{v + w + u + s) + \/\{u — u){v — v)\[Zi sgn(w — v) + Z~^ sgn{u — u)] 



= -(w + w + M + s) + sgn(w - v)^J\{u- u){v - v)\{Zi + Z2) 

= -{v + w + 11 + s) + sgn(u - u)yj\{u - u){v - v)\{Z^^ + Z^^). (41) 

In the thermodynamic limit L 00, all unimodular complex numbers are 
solutions for Zi . The maximum of the real part of E depends on whether all of 
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the values of Zi are (unimodular) or not. In the first case, the relaxation time 
is independent of the injection and extraction rates at the boundaries. This is 
called the ^ast phase. In the second case, the relaxation time does depend on 
the injection and extraction rates at the boundaries. This is called the slow 
phase. The terms fast and slow come from the fact that for fixed bulk reaction 
rates, the relaxation time in the fast phase is smaller than that of the slow 
phase. Now the physical parameter space can be divided to three regions, with 
different behaviors for the relaxation-time: 

a + a! < A, a + a' - (6 + 5') < ^ - region I 
b + b' <B, a + a' - {b + b') > A - B, region II (42) 
otherwise, region III 

where 

A u + s — \/\(u — u){v — v)\ — u + s ~ \/\{u — u){v — v)\, 

B := V + w — y/\{u — u){v — v)\ = V + w — \/\{u — u){v — v)\. (43) 

The relaxation time t is then 

{[v + w + a + a' + {u — u){v — v){a + a' — u — s)~^]~^, region I 
[u + s + b + b' + {u — u){v — v){b + b' ~ V — w)^^]^^, region II 
{v + w + u + s — 2Ke[^/JJ^^^^)(v^^^]}~^ , region III 

(44) 

It is seen that in region I, the relaxation time depends on the injection and 
extraction rates at the beginning of the lattice; in region II it depends on the 
injection and extraction rates at the end of the lattice; and in region III it 
depends on none of the injection or extraction rates. 
The details of the calculations can be seen in j2Dl ■ 



4 An asymmetric generalization of the Glauber 
model at zero temperature, as an example 

The Glauber model jSH] is a model for the relaxation of an Ising model towards 
its equilibrium with a heat bath at temperature T, which is based on the prin- 
ciple of detailed balance. This model is based on a three-site interaction: a spin 
is flipped with the rate /x := 1 — tanh /? J if the spins of its neighboring sites are 
the same as itself; it flips with the rate A := 1 + tanh/3 J, if the the spins of the 
neighboring sites are opposite to it: and flips with the rate 1, if the spin of the 
neighboring sites are opposite to each other. Here (3 = (fee T)~^, and time has 
been rescaled so that one of the rates is 1. 

It has been shown in |34| that at zero temperature, the Glauber dynamics 
is effectively a two-site interaction, in which two opposite neighboring spins 
become the same, with the rate 1 (independent of which is up and which is 
down, and which one flips). 
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An extension of the Glauber dynamics at zero temperature is a system with 
the following dynamics. 

with the rate u 
with the rate v 

with the rate v , , 

(45) 

with the rate u 

where A denotes a spin up (or a particle) and denotes a spin down (or a 
vacancy). In the ordinary glauber model, the system is left-right symmetric 
and u — v. Now consider a system with the dynamics (|45ll . on an L-site lattice. 
Assume also that there are injection and extraction rates at the first and last 
sites, with the rates a, a', 6, and h' introduced in the previous section. It is 
easily seen that for this system, of the eight bulk-rate parameters introduced in 
the previous section, only u, w, s, and w are nonzero, and 




From this, it is seen that 



w — V. (46) 



V 

X 



U + V 

y-^- (47) 

U + V 

So X and y are nonnegative and x + y = 1 . This means that there is a static 
phase transition, which occurs at u = v, the point corresponding to the ordinary 
Glauber model. 

For the dynamical phase transition, it is seen that A and B introduced in 
(O, are 

A := u — \/uv, 

B:=v- V^. (48) 

Of A and B, only one can be positive, and this is when u ^ v. So, as the 
injection and extraction rates are nonnegative, at most two of the regions I, II, 
and III may exist: for u > only the regions I and III; for u < v, only the 
regions II and III. The expression for the relaxation time is 

{[v + a + a' + u v{a + a' — u)^^]^^, region I 
[u + b + b' + uv{b + b' - vy^]^'^, region II (49) 
{v + u — 2^u v)^^, region III 

As a final note, it is easily seen that changing the dynamics H45|l into the 
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following 




0A, with the rate A 
AA^ with the rate u — X 
00, with the rate v — X 



A0, with the rate A' 
AA, with the rate v — X' 
00, with the rate u- X' 



(50) 



which contains diffusion as well, does not change the evolution equation of the 
one-point function. So the results of this section are valid for this interaction 
as well. The only difference is that the introduction of diffusion, restricts the 
parameter space of u and v, since the rates should be nonncgative. So u,?; > 



5 Extension to longer- range interactions: dy- 
namical phase transition of kinetic Ising model 
with boundaries 

It was pointed out in the previous section, that the Glauber dynamics at nonzero 
temperatures consists of a three-site interaction. Let us rewrite the reactions: 



where an upward arrow means spin up, a downward arrow means spin down, 
and A and fi are defined through 



as in the previous section. As it is noted in |33| . the evolution equation for the 
expectation value of the spins (each spin is either 1 (upward) or —1 (downward)) 
is closed. This is of course through for the spins at the bulk of the lattice. At 
the boundaries, one should introduce some two-site interactions and write the 
write separate equations for (si) and (sl). The new interactions introduced at 
the boundaries are 



A, A'. 



T T T ^ T i T and i i i ^ i t i 
T i T ^ T T T and i T i ^ i i i 
T T i ^ T i i and i i T ^ i T T 



with the rate /i, 
with the rate A, 
with the rate 1, 



(51) 



A := 1 + tanh/3 J, 
/X 1 — tanh (3 J, 



(52) 



T i i with the rate 5i , 

t I T with the rate 172, 

I T T with the rate 33, 

i T i with the rate 34, 



(53) 
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For the spin flip of the first site, and 

i i i with the rate hi, 

t T i with the rate /i2, 

t t T with the rate h^, 

i i-^ i T with the rate h^, (54) 

For the spin flip of the last site. It can be seen that the evolution equation for 
the one-point functions is closed, iff 

9i + .94 52+53, 

hi + hi^h2 + h3. (55) 

Provided this is through, the evolution equation is 

{si) = - 2{si) + {{si+i) + (s,_i>) tanh(2/3 J), l<i<L 
(si) = - (52 + 53)(si) + (51 - .92)(s2) + (53 - 51), 

(sl) - - (/i2 + h3){sL) + [hi - h2){sL-l) + (/l3 " h^). (56) 
The stationary solution for this is 

{sj) = 01x1 + 02 z^-^-\ (57) 

where 

zi = za"^ = tanh(/? J). (58) 

Oi and O2 are easily obtained at the thermodynamic limit L ^ 00; and it is 
seen that they are both smooth functions of the rates. So there is no static 
phase transition. 

For the dynamical phase transition, similar to the method used in the previ- 
ous sections, one writes the eigenvalue problem for the homogeneous equation. 
That equation reads 

E Xj = -2Xj + tanh(2/3J)(Xj+i + j 7^ 1, L, 

EXi^ -{g2 + 53)^1 + (51 - 52)^2, 

E XL^-{h2 + h3)XL + {hi-h2)XL-i. (59) 
The solution to this is an expression like (j^Hl, with 

S = -2 + tanh(2/3 J)(2j + 2~^). (60) 

and 

zl~^[2 - 52 - 53 + z{gi - 52 - tanh(2/3 J)) - tanh(2/3 J)] 
X [2 - /i2 - hs + z,{hi - /i2 - tanh(2/? J)) - z^^ tanh(2/3 J)] 
-zl-\2 - 52 - 53 + z-^{gi ~g2- tanh(2/3 J)) - z, tanh(2/3 J)] 
X [2 - /i2 - /13 + Z;"\ft.i -h2- tanh(2/? J)) - z^ tanh(2/3 J)] 0. (61) 
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Again in the thermodynamic Hmit L ^ oo, if all of the solutions to the above 
equation are unimodular, the relaxation time would be 

T = [-2 + 2tanh(2/? J)]"\ (62) 

which is independent of the interaction rates at the boundaries. If there are 
some solutions that arc nonunimodolar, and among them there is a solution 
with Ke{zi + z~^) > 2 (for the ferromagnetic case J > 0), then the relaxation 
time becomes larger than H62|l . This is the slow phase. The criterion for this, is 
that one of the roots of Ht)l|l passes the boundary 

y = ±i.-i)^^, (63) 

where x and y are the real and imaginary parts of that root, respectively. 

In the thermodynamic limit L — > oo, the equation (|61|l for roots with mod- 
ulus greater than one simplifies into 

{2 - 32 - 93+z^^[9i - 52 - tanh(2/3 J)] - z tanh(2/3 J)} 
x{2- /i2 - h3+z~^[hi -h2~ tanh(2/3 J)] - z tanh(2/3 J)} = 0. (64) 

So the criterion for the transition to occur, is that a complex number z — x + iy^ 
satisfies ()63|l and H64(l . The relation between rates, for this to occur, can be seen 
to be 

2[1 - tanh(2/? J)] - 92 ~ 9i = 0, or, 

[4tanh(2/3 J) - 2 + 52 + 53] [.9i ^ 52 - tanh(2/3 J)] + (2 - 32 - 53) tanh(2/3 J) = 0. 

(65) 

This is for the case when the first factor in (|64(l vanishes. A similar criterion 
come from the vanishing of the second factor, with the roles of giS and hi& 
interchanged. 

The detailed calculations can be found in |32| . 

6 Extension to multi-species systems: phase tran- 
sition in an asymmetric generalization of the 
zero-temperature Potts model 

A simple extension of the asymmetric generalization of the Glauber model at 
zero temperature, in troduced in section 4, is an asymmetric generalization of 
the zero-temperature Potts model. In this system, any site can have p -\- \ 
states (rather than the two states of the Glauber model). The bulk reactions 
are written as 

J AaAa, with the rate u 
yApAfj, with the rate f 
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Compare this with H45|l . One adds to this, reaction rates at the boundaries: 

A/j — > Aa with the rate A"^ at the first site, 

A/j Aa with the rate V^p at the last site. (66) 

For a (3, A"^ and F"^ are rates, and should be nonnegative. The diagonal 
elements of A and F are chosen so that 

sA = sF = 0. (67) 

Using these, the evolution equation for the one-point functions is written as 

(hj) = -(u + u)(nj) +u(nj_i) +v{nj+i), j ^1,L 
(hi) = A(ni) - w(ni) + w(n2), 
(hi) = F(nL) - u{nL) + u{nL-i). (68) 

6.1 The static phase transition 

The stationary solution to (|68|l can be written as 



(n,> = B - 
or 



(n,> = B 



«Qy-\ (ro, 



where 



sB = 1, 

sC' = sC" = 0, (71) 

and B and C (or C") must also satisfy conditions coming from the evolution 
equation 168|l at the boundaries. For m < w, it is better to work with ()69|l . Then 
it is seen that in the thermodynamic limit L oo, one has 

FB = 0, 

(A- w + m)C' = -AB. (72) 

The first equation has at least one non-zero solution for B, as 7 does have a left 
eigenvector s, with the eigenvalue zero. Depending on the degeneracy of this 
zero eigenvalue, there are one or several solutions for B. One should then put B 
in the second equation to find C The conditions on A, ensure that the real part 
of the eigenvalues of A are nonpositive, so v — u cannot be an eigenvalue of A. 
Hence, the second equation has one and only one solution for C (corresponding 
to each solution for B). 

For u > V, one uses H7(Jf) and (in the thermodynamic limit) arrives at 

AB = 0, 

{T -u + v)C" = -TB. (73) 
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Similar arguments hold for the solution of this equation, with the roles of A 
and r, u and v, and C and C" interchanged. Now apart from the question 
of the uniqueness of the stationary profile (which depends on the degeneracy 
of the matrices T or A) one can see that if u < u, the profile is flat near the 
last site, while for u > w it is flat near the first sites. This is the static phase 
transition, similar to what observed for the asymmetric generalization of the 
zero-temperature Glauber model. 



6.2 The dynamical phase transition 

The eigenvalue equation for the homogeneous part of (|68|l is 



(m + t;)Xj + mXa;_i + wXfe+i, j^l,L 



E Xi =A Xi — V Xi + V X2, 
-EXl =rXL - uXl + uXl_i, 

the solution to which is 

Xj = b + c Z2 , 

where ZiS satisfy 

E = ~{u + v) + V Zi + u . 
The second and third equation of (|74|l . can be witten as 

-{u + K)Z + ^Juv -{u + K)Z^'^ + ^/av 



-{v + T)Z^ + - 
where the variable 



>Z-^+i -{v + T)Z 



-L 



/uv Z 



-L-1 




Z, 



has been used, in terms of which 

E = -{u + v) + Vu^{Z, + Zri). 



(74) 

(75) 
(76) 

0, (77) 

(78) 
(79) 



H77|l has nontrivial solutions for b and c, iff the determinant of the matrix of 
the coefficients vanishes: 



'U V 



det 



-{u + A)Z- 
-{v + T)Z^ + ^/^Z^+^ 



-{v + T)Z-^ + ^JWliZ-^-^ 



(80) 



Again, if in the thermodynamic limit L ^ 00, for all of the solutions of the 
above equation Kc(Z + Z^^) < 2, the relaxation time is independent of the 
reaction rates at boundaries, and the system is in the fast phase. In order that 
the system be in the slow phase, it is necessary (but not sufficient) that at least 
one of the roots of H80|) has modulus larger than one. For this solution, H80(l is 
simplified in the thermodynamic limit to 



det[y/Tl^ Z-{v + r)] det[VT!^ Z - {u + A)] = 0. 



(81) 
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So, denoting the eigenvalues of A and F with A and 7 respectively, This root 
must be 

--" + ^ or (82) 



The real parts of the eigenvalues A and 7 are nonpositive. It is then seen that 
if M > w, then the real part of Z in the first case is not greater than 1, and 
Re(Z + Z^^) < 2. So in this case, the only relevant equation for finding the 
system in the slow phase is the Z = {u + X)/^/uv. A similar argument shows 
that ii V > u, then the first equation of 1)82(1 is relevant. So, without loss of 
generality, let's take u < v. In this case, b should be an eigenvector of T. 
However, one also demands 

sb = 0. (83) 

If b is the only eigenvector of F, with zero eigenvalue, then H83() cannot be 
satisfied, as s would also the only left eigenvector of F with zero eigenvalue. So, 
from the eigenvalues of F, one should put aside one zero eigenvalue, and consider 
only the other ones. Of course if the zero eigenvalue is degenerate, then 7 can 
still be zero. 

The system undergoes a dynamical phase transition, when x and y (the real 
and imaginary parts of Z) pass the curve 



y = ±(x-l)^^, x>l. (84) 

In terms of the eigenvalues of F, it is seen that the system is in the slow phase 
iff 



|Im(A)| < [Re{\) + u - \ hrr^^'^riT^. ' or Re(A) >2V^-w. 

y Zy/u V — Re(A) — u 

(85) 

A simple way to induce the phase transition is to multiply the matrix F by 
a parameter r. This means multiplying the rates of the reaction at the first site 
by r. As Re(7) < 0, one can see that for a large enough value of r, the value of 
Ke{'j)+u — y/uv will be negative (provided Re(7) ^ 0, that is, provided the zero 
eigenvalue of the matrix F is not degenerate) . So the system will be in the fast 
phase. It is also seen that as r tends to zero, either — Re(7) — u becomes 

negative, or in the first inequality in (|85|l the right-hand becomes greater than 
the left-hand side (which tends to zero). So, the system will be in the slow phase. 
Roughly speaking, increasing the reaction rates at the boundaries, brings the 
system from the slow phase (relaxation time depending on the reaction rates at 
boundaries) to the fast phase (relaxation time independent of the reaction rates 
at boundaries). A similar argument holds, of course, for the case u > v and the 
eigenvalues of the matrix A. 

The details of the calculations can be found in |31|. 
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